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A model is developed for the hadronic interaction in the two-nucleon system above pion threshold 
which is based on meson, nucleon and A degrees of freedom and which includes full meson retardation 
in the exchange operators. For technical reasons, the model allows maximal one meson to be present 
explicitly. Thus the Hilbert space contains besides NN and iVA also configurations consisting of 
two nucleons and one meson. For this reason, only two- and three-body unitarity is obeyed, and the 
model is suited for reactions in the two nucleon sector, where only one pion is produced or absorbed. 
Starting from a realistic pure nucleonic retarded potential, which had to be renormalized because 
of the additional 7r and A degrees of freedom, a reasonable fit to experimental NN scattering data 
could be achieved. 
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At present a very interesting topic in the field of medium energy physics is devoted to the role of effective degrees 
of freedom (d.o.f.) in hadronic systems in terms of nucleon, meson and isobar d.o.f. and their connection to the 
' underlying quark-gluon dynamics of QCD. For the study of this basic question, the two-nucleon system provides an 
— i , important test laboratory, because it is obviously the simplest nuclear system for the study of the nucleon-nucleon 
interaction in NN scattering including deuteron properties, and, furthermore, for testing this effective description in 
other reactions on the deuteron, for example, in elastic and inelastic electron scattering, in photodisintegration, and 
meson photo- and electroproduction. Moreover, due to the lack of free neutron targets, reactions on the deuteron are 
also an important tool to test our present understanding of neutron properties. As an example for the latter aspect, 
we would like to mention the recently much discussed determination of the electric form factor Ge,u of the neutron 
in d(e,e'n)p and d(e,e'ri)p or the investigation of the Gerasimov-Drell-Hearn sum rule for the deuteron ||. 
However, in order to separate the unwanted binding effects from the neutron properties, a precise knowledge of the 
O ' structure dependent effects is needed. 

Intense efforts over many past decades in experiment and theory have shown that for energies below pion threshold 
a satisfactory description of NN scattering and deuteron properties [|,||], as well as photo- and electrodisintegration 
, of the deuteron is achieved within the conventional framework of nucleon, meson and isobar d.o.f. although the 

description of observables is not perfect because for certain observables significant discrepancies remain unresolved 
as, for example, in elastic electron deuteron scattering ||. At higher energies, above pion threshold our theoretical 
understanding of the various experimental data is much less well settled. Even for deuteron photodisintegration (for 
a detailed review see [0 ) , none of the various theoretical approaches like the diagrammatic method of Laget 0, the 
framework of nuclear isobar configurations in the impulse approximation [fiof , the unitary three-body model [fllf or 
the coupled channel approach (CC) jl2],[l3| is able to describe in a satisfactory manner the whole set of experimental 
data on differential cross sections and polarization observables for energies covering the whole A resonance region. 

A common feature of most of these approaches is the extensive use of the static limit for the meson propagator 
which enters the hadronic interaction and the electromagnetic two-body exchange current operators, although there 
is little justification for that in view of the relatively high excitation energies involved. Also from the point of view 
of special relativity, a non-instantaneous interaction would be required. The reason, why this static approximation 
is still used, is the enormous simplification of the operator structure becoming local, and thus is much simpler to 
evaluate numerically. It has already been conjectured in ]l3| , that one main reason for the above mentioned failure 
of the theoretical description lies in the neglect of meson retardation in the meson exchange operators. Indeed, 
first results 0,^5| based on the thesis of M. Schwamb |l6| have shown that a much better description of deuteron 
photodisintegration compared to the static approaches is obtained if meson retardation is retained. 

In this paper, we want to present a realistic, but still tractable model for a retarded hadronic interaction within a 
nonrelativistic framework which is suitable to be used as input in electromagnetic reactions on the deuteron like, for 
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example, photo- and electrodisintegration or pion production. There exists already a variety of models for treating 
retardation in NN scattering based on three-body theory, for example, the work of Kloet and Silbar and Tanabe 
and Ohta In these calculations, however, the driving force is basically one-pion-exchange so that the lower partial 
waves in NN scattering cannot be described reasonably well, even at energies below pion threshold. The neglected 
short-range interaction is, however, of crucial importance for electromagnetic and hadronic breakup reactions on the 
deuteron. Tanabe and Ohta, for example, have "patched up" this problem in photodisintegration of the deuteron |Tl| l 
by mixing retarded and static frameworks, namely the static Paris potential for the final state interaction with respect 
to the nucleon one-body current and the it MEC, whereas for the A contributions they employ a three-body model 
in which heavy meson exchange has been neglected. Kloet and Silbar on the other hand have taken a retarded pion 
exchange and static heavy meson exchange in an improved model of NN scattering p9| . Also other treatments 
of retardation in the AA-interaction like, e.g. Refs. |2^,|23| have some drawbacks as will be discussed later. The 
shortcomings of these approaches have been overcome to a large extent in the present model which ensures, due to 
the use of the Elster- Potential (f(],|2l[] as an input, a reasonable description of all relevant partial waves below and 
above pion threshold within a consistent framework in time-ordered perturbation theory. 

Thus we will present here the basic framework of a model for the hadronic interaction in which retardation in the 
exchange operators is fully included. We would like to emphasize the fact that this formalism can be applied to any 
reaction on the two-nucleon system for excitation energies up to about 500 MeV in which not more than one pion 
is produced or absorbed. But in this paper we will restrict ourselves to the hadronic interaction in studying NN 
scattering and deuteron properties only in order to fix all free hadronic parameters of our model. The application of 
this formalism to other hadronic reactions like ird scattering as well as to electromagnetic reactions will be deferred 
to forthcoming papers. 

The conceptual basis and the main features of our model are laid out in Sect. [n]. In view of the fact that conventional 
retarded interactions in pure nucleonic space are energy dependent and thus non-hermitean, we enlarge the Hilbert 
space by considering explicitly meson and A d.o.f. in order to start with a hermitean hamiltonian. The retarded 
interaction is then generated by meson-nucleon and 7rAA vertices. For reasons of simplicity, we restrict ourselves 
in this work to configurations with only one meson present besides the baryons. Special attention is laid on the 
question of nucleon dressing and the corresponding renormalization of operators as well as the fulfilment of two- and 
three-body unitarity. With respect to the latter, we have incorporated the ird channel in a realistic model showing 
some interesting effects. The NN scattering T matrix is derived, and the structure of the deuteron is discussed. A 
field theoretical realization in the form of a one-boson-exchange model is developed in Sect. Ill with inclusion of the 
A isobar which is conceptually similar to the work of the Bonn group p3 but differs in some essential details. Then 
we will present and discuss in Sect. IV the results for AW scattering comparing them with experiment as well as with 
other theoretical approaches. Finally, Sect. ^ contains a summary and an outlook. 



II. BASIC CONSIDERATIONS 



A. The Hilbert space 



Our model for the description of a nuclear system allows besides n nucleons configurations with one additional 
meson present or where one nucleon is replaced by a A isobar. It is similar (but not identical) to the approach of 
Sauer and collaborators p3| , p5| -[2^] which is also used in [0,p8p3j. Thus the model Hilbert space is subdivided 
into three orthogonal spaces according to the different configurations containing either n bare nucleons (7i}^), n — 1 
nucleons and one A (Ti^), or n nucleons and one meson (7ij^'), i.e. 

«M = ftg] 0W N eW W. (1) 

The "bar" indicates a bare nucleon to be distinguished from the corresponding physical nucleon, denoted without 
a bar. This distinction is necessary in order to take into account the dressing of a bare nucleon by meson-nucleon 
loops. For the A isobar we will not make this distinction, because the self energy contributions from ttN loops will be 
retained explicitly, whereas the dressing of the bare nucleon to become a physical nucleon will be incorporated into 
the effective operators by dressing factors. In Ti}^ only one meson is present besides n nucleons, i.e., no components 
with two or more mesons are taken into account (one- meson-approximation). This limitation creates some pathologies 
as will be pointed out later. In detail it means 

W W= (2) 
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considering as mesons 7r, 77, a, 5, lo, and p. 

In order to distinguish the various sectors of H)- n \ we introduce corresponding projection operators by 

p N n [n] = n [ ™ ] , p A n [n] = n [ 2 ] , p x n [n] = , and p = p R + p A . (3) 

The latter is introduced since later we also will need the projection operator on the combined pure baryon subspace 
© H^. Moreover, in view of the different subspaces in the meson-nucleon sector, it is useful to decompose the 
projection operator Px into a sum of six orthogonal projectors corresponding to the six mesons considered 

Px = P * with P * Hln] = n * ] for x e ^' 1> CT < w ' p} • ( 4 ) 

x£{tT, 7], <T, 5, LO, p} 

Using the notation 

Q af} = PMPp , a,p£ {N, A, X} , (5) 
any operator f2 acting in Ti\ n ^ can be written as a symbolic 3x3 matrix 

"iViV ^NA ^NX 

= ( ^an ^aa Max | • (6) 

^XN ^XA &XX 



B. The Hamiltonian 



The Hamilton operator H of the model can be divided into a diagonal kinetic part and an interaction describing 
the emission and absorption of a meson by a baryon and, in addition, an ab initio baryon-baryon interaction 

H = Ho + Hi , (7) 

where the bar indicates that it refers to bare baryons. In particular, the nucleon kinetic energies are determined 
by the bare nucleon mass M N . However, due to the truncation of the Hilbert space with respect to the number of 
mesons, no explicit dressing of the bare nucleons is possible in Ti)^ © ■ This is one of the pathologies which one 
encounters in the one-meson-approximation. Therefore, we will use the physical instead of the bare nucleon mass in 

[nl [nl 

Tt A ©Hy, so that the formally suppressed meson-nucleon loops can be taken into account at least effectively by the 

physical nucleon mass. For simplicity, we use in Ti}^ the nonrelativistic expressions for the kinetic energies of nucleon 
and A isobar because we will treat the A nonrelativistically. This approximation is, however, not crucial. For later 
purposes we will introduce as kinetic energy in addition a second type of one-body operator Hq which differs from Hq 

[nl 

in the pure nucleonic sector only, referring to physical nucleon kinetic energies also in TL^ . Thus we have in detail 

n 

H 0,NN = (8) 

n 

H 0,NN = X) ftjV (•?')' ( 9 ) 

i=i 

n-1 

ffo,AA-ffo,AA = ^^(j) + ftA, (10) 



Ho, xx — H xx = J2hN(j) + h x , (11) 
with the various kinetic energies (with masses m x , x e {it, r], a, S, u, p}, M A , Mfj and Mn, respectively) 
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(N(p')\h N \N(p))=e N (p)S(p 
(N(p')\h N \N(p))=e N (p)5(p 
(N(p')\h n N r \N(p))=e% r (p)5(p 

{A(p')\h A \A(p)) = e A (p)S(p 
{x(q')\h x \x(q)) = (2ir) 3 2u x (q)5(q 



P) 
p) 
p) 

P) 



zn{p) 



M%+p\ 



M 2 N + P 2 , 



e n N r (p) = M N + 
e A (p)=M° A + 



P~ 



2M N 

P 2 
2M A ' 



, (p) = \/ ml + p 2 , for x e {it, r), a, 5, uj, p} 



(12) 
(13) 
(14) 

(15) 
(16) 



In view of the two choices for the kinetic energy, one obtains besides the interaction Hi, defined in (^) another 
interaction operator V° as defined via 

H = H Q + V°. (17) 
The various components of the interaction V° are depicted in Fig. [I]. First of all, V° contains, because of 

V° = H I + H -H , (18) 
besides the interaction Hi, a diagonal counter term 

V® =H -H , (19) 
which is nonzero in Tij™' only, and which is a pure one-body operator 



V tiN = V MvC?) ' With V tiN = h-N-h N . 



(20) 



3=1 



As already mentioned, we allow in addition in Vp P a two-body part Vpj^ which describes an ab initio given hermitean 
interaction between two baryons, i.e., 



V° 



M [2] 

/ pp , 



' pp — " pp 

which will be specified later. The nondiagonal components V XN and V XA are one-body operators 



(21) 



3=1 3=1 



(22) 



describing the emission of a meson by a baryon. The remaining interaction V xx consists in the one-meson- 
approximation in principle of two parts 



yO _ V 0N , \p yOx 

a;£{7r, rj 1 a, 6, uj, p} 



(23) 



The first one describes a meson-nucleon interaction with the other nucleons as spectators and the second one two 
interacting nucleons with a meson x and the remaining nucleons as spectators (see Fig. Q). In consequence, V xx is a 
one-body operator 



ON 
XX 



(24) 



In the previous work of Wilhelm et al. gHU^ and Bulla et al. ||, V xx was neglected completely for practical 
reasons. Above the %d threshold, however, this approximation leads to several problems, in particular for three-body 
unitarity (see the discussion in Sect. HID), which are solely due to the neglect of V xx . We therefore set 



and for x ^ n 



whereas only V xx will be retained nonzero. 



V 0N _ n 

v xx — u ) 



v xx = o, 



(25) 
(26) 
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C. NN scattering 



Now we will consider NN scattering in the two-nucleon sector. The renormalization of a free nucleon state \N; p) is 
sketched briefly in Appendix A. The scattering states of two physical nucleons | NN; p, a)^ in the cm. system with 
the asymptotic free relative momentum p, energy E^ N = 2cn{p) and a complete set of internal quantum numbers a 
are given by 

\NN; p, fl )( ± ) = N- 2] \p ) (l + G (E?" ± ie)T°(E» N ± iej) \NN; p, a) , (27) 

where Nn] (p) is a renormalization constant which appears because both, the bare as well as the physical free two- 
nucleon states are normalized to the 5 function. The transition amplitude T° satisfies the Lippmann-Schwinger 
equation 

T°(z) = V° + V°G (z)T°(z) (28) 

with the free propagator Go 

G (z) = (z - Ho)' 1 , (29) 

which is represented in Fig. |2[ As is shown in detail in the Appendix B, one finds for the T matrix element 

(NN; f, a'\T(E» N + ie)\NN; p, a) = N^(p)(NN; f, a'\T°(E^ N + ie)\NN; p, a) 

= (NN; f, a'\T^(E^ N + ie))\NN; p, a) . (30) 

Here, T con - the superscript "con" refers to connected diagrams (see Appendix B) - obeys a Lippmann-Schwinger 
equation 

T c P °P(z) = R(z) Vgp{z) R(z) + R{z) Vfi?(z) R{z) G Q (z) Tg?(z) . (31) 
Its driving term contains a "renormalized" interaction 

V&(z) = (Z^y 1 V° P con (z) (Zffi)- 1 , (32) 

where Vpp° n , as defined in Appendix B, comprises besides irN loop contributions to the A self energy the genuine 
retarded baryon-baryon interaction. Its diagrammatic representation is shown in Fig. ||. Furthermore, the "dressing 
operator" is given by 

R(z) = Z$Z^(z). (33) 
Here the two-body renormalization operator Zy 2 \ is defined as 

Zfa(z) = l + J dz'S(z'-H ) [V% x G (z')G (z)V^] d . g , (34) 

where the subscript "dis" refers to disconnected diagrams (see Appendix B), and which differs from unity in WW 
only. Its onshell value is 

(%) a = 1 + f dz'S(z' - H ) [V% x G (z') Go(z') V xS ] dis . (35) 



Thus NN scattering is unambiguously fixed by the onshell matrix element of TS 

For later applications we need also the offshell form of Tp°p(z) for which one finds the following expression 



^con 
~ NN' 



Tp^(z) = Gq -\z) {G A) (z) + G ( A \z)T P p>(z)G ( Q A \z)} G \z) G^(z) . (36) 

The auxiliary amplitude Tp°p(z) is given by the integral equation 

T PP \z) = V$ PP (z) + V^p P (z)G { A \z)f^(z) , (37) 
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where the driving terms V^° PP (z) (i = 1, 2) are defined by separating in V P ° P (z) the irN loop contributions to the A 
self energy, subdividing it into two parts 

Vpp( z ) = Vfi PP (z) + V$ PP (z) , (38) 

where 

Vfi n PP (z) = [V AX G (z)V XA ] dis , (39) 

V$ n PP (z) = V P 2] P + [V PX G (z)V XP ] con + V PX G a (z)T x (z)G (z)V XP , (40) 

with the renormalized interactions 

^ = (%]r 1 ^l I (%])- 1 , (4i) 

V XP = (Z^)- 1 V XP (Zffi)- 1 . (42) 

Here, T x {z) is the NN scattering matrix in the presence of a spectator meson fulfilling 

T X (z) = V xx + V xx Go(z) T x (z) . (43) 

Its diagrammatic representation is shown in Fig. ||. Thus the term V^° pp (z) contains solely the intermediate irN 
loop contributions to the A propagator. Furthermore, we have introduced in ( |36| ) a "dressed" propagator 

G [A \z) = G (z) + G a {z)V^ PP {z)G {A \z) , (44) 

which takes into account the dressing of the A in TC^. In TC^ and H x \ Gq A \z) is identical to the free propagator, 

wher eas in H A one gets a simple connection between Gq A \z) and the A propagator g A (z) in the one- A sector (see 
Sect. HH) 

?( A ) f„\\MA. *\ - /A. ^/U. (~ 1\J.. 



(NA; p'\G$ A (z)\NA; p) = (A; p'\g A z - M N - ^— |A; p) , (45) 



where p is the relative momentum of the NA system, and its reduced mass is denoted by 

M° A M N 

MWA = Mi + M N ■ (46) 
In addition, we need the baryonic and mesonic components of the NN scattering states for which one finds 

P N \NN; p, a) (±) = 5^ (l + G (z)Tffi(zj) \NN; p, a) , (47) 

Z [2] 

P A \NN; p, a) {±) = G (A) (z)f%$(z)\NN; p, a) , (48) 
P X \NN; p, a) {±) = G x (z)(v XN Zffi + V XA ) \NN; p, a) , (49) 

where z = E^ N ± ie, and G x (z) describes the propagation of two interacting nucleons in the presence of a spectator 
meson 

G x (z) = (z-H . xx -Vxx)- 1 

= Ga(z) + Go(z)T x (z)G (z). (50) 



G 



D. The deuteron 



Due to its vanishing isospin, the deuteron cannot contain a NA component. Thus, the deuteron state \d) can be 
separated into a nucleonic and a mesonic component 

\d) = Px\d) + P x \d) , (51) 

which can be determined by the Schrodinger equation in the cm. frame 

(H + V°) | J) = M d \d) (52) 

with M d = 2Mn — £b as deuteron mass and Sb its binding energy. Eliminating the mesonic component, we introduce 
the purely nucleonic part Pft\d) as an effective renormalized deuteron state by 

\d) = Z [2] {M d )P n \d) , (53) 

with the normalization (d\d) = 1. It is straightforward to show that \d) obeys the equation 

(h + R(M d ) (V|l + [% G oWF^] OT ) R(M d )) \d) = M d \d) , (54) 

which contains only renormalized quantities. From the effective state \d) the original deuteron state \d) is obtained 
by 

\d) = i- ( + G (M d )V xN R(M d )) \d) , (55) 



N d \ Z° 2] 



with the renormalization constant Nd 



N 2 d 



(d\ (l-§- z {R(z) [Vf} X G (z)V x ^] con R(z)}\ z=M ) \d) . (56) 



Due to the absence of NA components and the vanishing of V xx in (t = 0) channels, the quantity VkL in ( p34] ) 
can be put equal to zero in retarded calculations. On the other hand, in static approaches, because of the choice 
V^ x = 0, one has to identify V^ffl with the chosen realistic NN potential V^ e § 1 (see the next section). 

III. FIELD-THEORETICAL REALIZATION 

Now, we will introduce a field-theoretical realization of the hadronic interaction. In view of various approaches in 
the literature, it will be useful to distinguish two types of realizations which differ in the treatment of the interaction 
in Hjj only, i.e., in V^^^ of JIo|) and which we will coin "static" and "retarded" approaches. In order to illustrate 

the essential differences, we will set V xx equal to zero for the moment being for simplicity. Then the pure nucleonic 
component V^Nn cons i s ts of two contributions: (i) a given hermitean, energy independent potential generated from 

, and (ii) a retarded one- meson exchange potential 

Vffi(z) = [V Rx G Q (z)V XN ] con . (57) 

In the following, we will use the notion static approach for the case that any explicit meson-nucleon vertex 
vanishes. Thus in the static case the retarded one-meson exchange vanishes identically and the NN interaction 

[21 

Vj^™ is generated by alone. One should note, however, that even in the static approach retardation is still 

contained in the interaction V£°\- Furthermore, there is no distinction between a bare and a physical nucleon, and 
the operators Zi 2 ] and R are both equal to the identity. Consequently, one can leave out the "bar" in the notation. 
Moreover, the upper index "0" in the interaction operator V° can be dropped. In order to have a realistic description, 
one then has to identify Vffl N = V^jy with a realistic NN potential model VJvIvm i- e -> 

V[2]nn( z ) = Vnn = ■ (58) 
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Furthermore, one has to keep in mind that due to the coupling to the NA and ttNN states, the realistic potential 
has t o be renormalized in order to avoid double counting of parts of the interaction as will be discussed below in 
Sect. [II E[ Such an approach has been used in |25|,^6| and also in the coupled channel calculation of |13 2{J. The 



obvious advantage of this framework is its simplicity. For excitation energies up to about 500 MeV, only pions could 
be created via the ttN A vertex. 

In the retarded approach, on the other hand, one chooses ^ and VSjS = 0. Neglecting V^ x and V?v> for a 
moment, the NN interaction is generated completely by the retarded one- meson exchange part V^ e ^(z). In order to 
obtain a realistic description, one then has to consider besides the pion also heavier mesons explicitly. Otherwise, the 
NN interaction would be generated by a retarded one-pion exchange potential only, which would result in a rather 
crude description of experimental data. In the present work, we use the potential of Elster et al. j2Cj]2l[] . Also in this 
case a renormalization of the TV TV potential will be needed if one includes explicitly the iVA and irNN channels. 

It is obvious that besides these two extremes various alternatives are possible. For example, Bulla and Sauer 0] 
used in their extension of the original model |25|,|26| the choice 

Vk , Vg x = Q , x G {p, u, a, 5, v } , v|L = V&f - V£% (z = 2M N ) , (59) 

and for the Paris potential PQ], which has been ren ormal ized with respect to n exchange in order to avoid 



double counting as will be discussed in detail below in Sect. HIE 



A. The component V^ x 

The components V? „ and V®^ = (V^) 1 ^ are explicitly present in retarded calculations only. For the meson- 
nucleon vertices, we have taken the usual couplings for pseudoscalar, scalar, and vector mesons whose explicit forms 
we have taken from p4[ ]. At each vertex we have furthermore introduced a phenomenological hadronic form factor 
F x (q 2 ) parametrized in the conventional monopole (n x = 1) or dipole (n x = 2) form 

/ A 2 _ 2\ ». 

F x {f)= [ A 2 +g ^ j ^ xe{n,r,,a,5,oj,p}, (60) 

where the cutoff parameters A x are treated as free parameters to be fixed by fitting the NN scattering data below ir 
threshold and the deuteron properties. 



B. The coupling V& x and the dressing of the A isobar 

In view of the strong coupling of the A isobar to the irN system, we restrict ourselves to the coupling of the A to 
the ttA^ channel, i.e., ^ only, for which we take the usual nonrelativistic form for the one-body vertices of ( p^ ) 

fO 

^{p'X'-JM^) = ^^^Xa.a'^ajv ■^Xn,\F^n-k{^ 2 ) t an^, (61) 

where \A, y an d X^v a denote the nonrelativistic A and nucleon spinors, respectively, and the spin and isospin 

transition operators c? A jv and r AJ v are fixed by the reduced matrix elements (§||^lr|||) = (f ll 7 /^!!^) = 2- Again a 
phenomenological form factor 

F ^?)=(% MJL ^Y* (62) 



has been introduced. Similar to |13] , f26[ , the free parameters f% Nlr and Aaatw are fixed by fitting nN scattering in the 
P33 channel. One finds for the dressed A propagator g&, depicted in Fig. B, 

g A (z) = (z-Ml -Ea(z))- 1 (63) 

with the A self energy 

Z A (z)=v° AlT go(z)v° A . (64) 
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Because of our choice v xx = and the one- meson-approximation, background mechanisms like, e.g., the Chew-Low 
term have to be neglected as in 0,p6|. In detail, one obtains for the A mass and width, respectively, 



4tt 

M A (W) = M° A + —P 



dq' q' 



./A 



ANt, 



(27r)32w„(g') \ m* J W - ^{q 1 ) - e%{q') 



T A (W) = < 67T u J „(q)+M N { m x J 

o 



f 1nM 2 ) fOT W>m n +M N , 
for W <m^+ M N , 



(65) 
(66) 



with the invariant energy W = ^^(q) + e r ^f(q). For the free parameters /Xtvtt' ^an-k and M A , our fit to the solution 
SM95 of M yields 



(/aJVtt) 



4tt 



0.9452 , A 



ANir 



482.11MeV, M A = 1281.7 MeV, n Av = 2 . 



whereas in ]l3| , p6t 



(/aATtt) 
47T 



= 1.393 , K AN k = 287.9 MeV , M A = 1315 MeV , n A ^ = l 



(67) 



(68) 



has been used. We would like to emphasize that the values in (67) have been obtained by fitting simultaneously ttN 

(3/2) 

scattering in the P33 channel and the MW multipole of photopionproduction. The reason for this procedure will 

become apparent in a forthcoming paper on e.m. reactions on the deuteron, because it turned out that a reasonable 

(3/2) 

description of the most important MW multipole in the A region is not possible in our approach if only irN 
scattering is considered for the fit of f AN ^, Aajvtt and M A . 



C. The interactions TVTV <-> TV A and TV A <-» TV A 

The time ordered diagrams of the interactions TVA <-> TV TV and TVA <-> TVA are depicted in Fig. EL In view of the 
truncation of the model Hilbert space, not allowing explicit meson- TVA or meson- AA configurations, the contributions 
to the three diagrams (d) through (f ) have to be described by the ab initio potentials V^W , VjL A and V AA in the 
energy-independent limit, represented by the diagrams (d') through (f) in Fig. |]. In the static approach, the A-TV 
mass difference is neglected in the meson propagator of the diagrams (d') through (f) while it is retained in the 
retarded approach. As in the work of |23] , |26|| , we consider ir and p exchange only. With respect to the other three 
time-ordered diagrams (a) through (c) of Fig. EJ, for which we consider only it exchange, one has to distinguish 
retarded and static approaches. Whereas retardation is kept for diagram (c) in both cases, diagrams (a) and (b) are 
only considered in the retarded framework. Otherwise they are contained in V^j^ and V^^. 

The corresponding matrix elements for retarded tt exchange have the following structure 

/tt/aWtt t NnO-) ' 7iVA(2) 



(TVTV; p'\[V£ x G (z) V£ A ] co „|TVA; p) = F„{q 2 ) F AN „{q 2 ) ■ 



ml (27r)32^(g) 



x fa(i)-?)(y-r> (1h2) (69) 



(TVA; pWZx Go(z) V XA ] con \NA; p) = F^q 2 ) H^'/*?^ 

„ (^Ajv(I) -<f) (%a( 2 ) • <f) 



e%(p)-ej(p>)-u n (q) 



(1 «-> 2) , (70) 



where q = p' — p. For the evaluation of these expressions, the nonrelativistic reduction of the vertices has been used. 
As in (|65|), we take the nonrelativistic nucleon energies e^{p) in the ttNN propagators for simplicity. 

The remaining one-pion exchange diagrams (d') through (f ') in Fig. [| and the energy independent limit of all p 
exchange diagrams corresponding to (a) through (f) in Fig. EJ are included in V^?', V A ft and V^j^', respectively, 
which we decompose into tt and p exchange parts according to 
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Vxl ] {z) = V°Z\z) + V^\z), (71) 
v2% ] (z) = V°£\z) + V^\z). (72) 

Explicitly, we take for them in the retarded approach the following expressions 

/ at at -*i I t ^0 (x) I at a ->\ p /-2\ s / ->2 \ fx /ajvit T JViv (1 ) ■ "Tzv A (2) (o'jvjv(I) ' ?) (<?/va(2) • ?) n /_„-. 

<M;P| ^ I^A;p> = F )r (« )2W« )— j (27r)32 ^ (g) Mj v~M A -^( g ) + C 1 ~ 2 ) ■ ( 73 ) 

p-|^|iVA; # = F^ N (q 2 ) ^ f^^ ^gf^W^rf + (1 - 2) . (74) 

mi {2TT) i 2LO^{q) 2M N - 2M A - uj n (q) 

The corresponding expressions for the static approach for are obtained by neglecting here the N-A mass 

difference. The form factor i^y/vV) parametrized as in ( |62] ) with parameters A^attt and n An> and the coupling constant 
f^Nir are n °t identical with F^n^ and f ANll , respectively, which were fixed by fitting irN scattering data, where irA 
configurations are not considered. Consequently, ^aattt and f A pf v can t» e treated in principle as free parameters to 
be fixed by NN scattering (see the following section). 

With respect to p exchange, we consider besides the usual pNA interaction density Sj 



(75) 



where ip A , tp, and </)„ denote the A Rarita-Schwingcr spinor, the nuclcon Dirac spinor, and the p meson field, 
respectively, an additional alternative 
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AM 2 N 



^J 5 r A ^-(d^)(dJ u -dJ^+h.c. , (76) 



which, to our knowledge, has not been discussed in the literature. A nonrelativistic reduction yields in the cm. frame 
(NN; ?\V$> \NA; p) - ' ^ a(2) ( 1 — 



XA ' (2nf2cj p (q) \M N -M A -uJ p (q) u p {q) 

o 
p 

2M N 



J A, 

m 

{ffjfifO) x q) ■ (ctjv A (2) x q) 



xF P (q>)F ANp (q*) \ -^4^ A (2) • (*x P) 



2M N 

a ANp (^(2) x q) • (^ A (2) X q ) 1 + (1 <-► 2) , (77) 



2Mn 

, 2 



(p) I A ,rv -,\ _ U%Np) T A ^(1) ■ T NA (2) 



(AN;p'\V^'\AN;p) 



AA I— >*/-— ^2 



p (27r) 3 2o; p ((7) \2M N -2M A -u p {q) u p (q) 

2\ 



with 



and 



xF 2 ANp (q 2 ) (? AJ v(l) x g-) • (<?tf A (2) x g) + (1 <- 2) (78) 



p'-p, P=l(p'+p) j (79) 



a A jv P = l--— . (80) 

4M iV / AWp 

Again the form factor F A M P (q 2 ) and coupling constant / AWp are fitted to NN scattering, whereas / p are g° p are fixed 

by the values of the Elster-potential. Note that the terms proportional to (q x P) and (<Tjvjv(2) x 9*) ■ (^jva(2) x q) 
are usually neglected in the literature, for example in 23 2q]. 
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Finally, with respect to the static approach, only the NA potential (|7(]) , represented by the diagram (c) in Fig. [|, 
is generated by iteration of the irNA vertex, as has already been mentioned above. All other diagrams of Fig. |] have 
to be described by Vpp^ where in addition the N-A mass difference is neglected. Consequently, we obtain in this 
limit 

(NN; p'K^\NA; p) = -F^)F AN ^) f M^ * N »^ ^ A . (2) 

mi (27r) 3 ^(<j) 

x(<? A , w (l) ■ q) ■ {a NA {2) ■ q) + (1 «-> 2) , (81) 

whereas Eqs. (Q), ([77|), and (JT^) apply also to the static limit except for the neglect of the N-A mass difference. 
Note however, that in this case all coupling constants and cutoffs in V XA and can be treated in principle as free 
parameters due to the choice V^ x = 0. 



D. The interaction V xx 

The simplest choice for the diagonal interaction in the subspace W-^ would certainly be V xx = 0. However, 
such a choice would lead to severe inconsistencies with respect to pion d.o.f., because it would lead to a violation of 
three-body unitarity. The reason for this violation lies in the fact that for excitation energies up to about 500 MeV, 
ttNN and nd states can exist in this sector as asymptotically free states (2ir and 3w states are not allowed due to the 
one-meson- approximation). It is therefore obvious, that the interaction V xx , which describes the NN interaction in 
the presence of a spectator pion, must be considered at least in the 3 Si- 3 -Di channel which we will henceforth refer 
to as the ird channel. Otherwise, the ird state would not be present formally and reactions like 7rd — > nd could not 
be studied without creating inconsistencies. In p7|-p9[, for example, reactions with a nd state in the initial and/or 
final state were studied without considering V xx ^ leading to a violation of unitarity above the ird threshold. On 
the other hand, a nd state could not be generated if V xx — 0. Therefore, V xx has to be nonzero, since we are 

[21 

interested in the construction of a unitary model up to the 2ir threshold. All other diagonal interactions in H x do 
not affect three-body unitarity for energies up to 500 MeV and, therefore, can safely be set equal to zero for the sake 
of simplicity (see (|2^) and (|2^)). We consider this choice as the minimal requirement in order to satisfy three-body 
unitarity. 

Thus we retain in V xx solely the interaction of two nucleons with isospin t — in the presence of a pion as spectator, 
where we restrict ourselves to an interaction, called V d , which acts only in the 3 <Si- 3 Di channel. Then V xx can be 
written in the form (see Fig. 



^ = / wftw) k( * )>nMd + ' (82) 

where |7r((f)) denotes the spectator pion state with momentum q. For practical reasons, we use for V d a separable 
interaction of rank 1 



which satisfies the Schrodinger equation for the nucleonic component of the deuteron 

(H + V d {M d ) - M d ) \d) = . (84) 
The separable structure of V d leads to a rather simple expression for the relevant amplitude T x of ([l3|). One obtains 

f q , ; h(q ))T d (M d + z - 

(27r) 3 2w 7r (q) 1 yq " V 4M N ' 4M N 

where the amplitude T d is given by the analytic expression 

T d (a r) = Gz\a)md\GZ\<*) (m 

V ' j (d\Go 1 (a)\d)-(d\Go 1 (a)G (z)Go 1 (a)\d)- 



The shifted arguments in T d in ( p5[ ) have their origin in the fact that the two nucleons, interacting via V d , have a 
total momentum —q. We use the nonrelativistic energy in order to separate the cm. energy exactly. It is therefore 
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natural to use the nonrelativistic e nerg y also in the propagator Go in ( pq ) . Boost contributions to the deuteron wave 
function are expected to be small [[32| and thus are neglected. In our numerical evaluation, the intermediate ttNN 
propagation, entering into the term [Vpx Gq(z) T x (z) Gq{z) Vxp] C on (cut A and B in Fig. |?]) of Eq. (EtJ) is treated 
in the static limit for simplicity. Moreover, we use for the nN vertex V^ x in [Vp X Gq(z) T x (z) Go(z)V X p] con the 
nonrelativistic version. Both approximations lead to a slight violation of unitarity which, however, is not critical. 

At the end of this subsection, we will discuss briefly the quality of the separable interaction in the 3 Si- 3 Di channel. 
For the deuteron pole, i.e., for W — Md, the resulting amplitude T d is identical to the exact solution obtained with the 
Bonn-OBEPR and Elster-potentials, respectively. However, for W ^ Md this equivalence breaks down. Concerning 
the resulting phase shifts, one obtains a strong deviation from the exact calculation for the 3 Di phase shift and the 
mixing parameter e (see Fig. ^). The most important 3 Si channel, however, is described reasonably well. These facts 
indicate the limits of the separable ansatz. 



E. Renormalization of the realistic NN potential 



A problem of double counting in the nucleon-nucleon interaction appears if one starts from a potential incorporating 
effectively certain d.o.f., which have been projected out before but which are introduced again explicitly later on. For 
example, let us consider a realistic NN potential like Bonn-OBEPR, Bonn-OBEPQ HQ, Paris Q, Argonne V u @, 
or Nijmegen [Q. These potentials act in pure nucleonic space and are fitted to deuteron properties and NN scattering 
data below ir threshold. However, if such potentials are used in a model with explicit A d.o.f. within a coupled channel 
approach, the problem of double counting becomes evident, because, for example, the dispersive box graphs depicted 
in Fig. ^ and implicitly present already in Vnn would be considered explicitly in addition. 

A simple way out of this problem is the box renormalization of Green and Sainio [ |35[ which consists in a subtraction 
of the box diagrams at a fixed energy Eq from V<r'^° n in (B7). We will adopt this method also here and use - similar 
to previous work |l3|,^3],^5 J26 2jJ - the value Eq = 2Mx. With respect to the above discussion we will distinguish 



(i) static calculations (VW X = 0): 



V, 



0[2] 



NN 



V'real 
NN 



V° A G (z)Vl N \ 



V° x Go(z)T x (z)Go(z)V XN \ 



(87) 



and 



(ii) retarded calculations (V$ ^ o, Vfi% 1 = 0) 



V N f=-V^ ll (z)G (z)V?( ull (z) 



z=E 



V£ x G (z)T x (z)Go(z)V x 



XN\ 



(88) 



where we have defined 



V°l uU (z) = V£ A (z)+[V? x Go(z)V : 



x a] con 



(89) 



It should be emphasized that the intermediate NA and Ttd states have isospin t = 1 so that the box subtraction 
has no influence on the (t = 0) channels, especially for the deuteron. However, such a subtraction would appear if 
also AA configurations would be included explicitly. Furthermore, we would like to remark, that this subtraction is 
not a fundamental ingredient of our approach. It is only a relatively simple recipe in order to incorporate a given 
realistic NN potential, which does not contain explicit A degrees of freedom, into a NA coupled channel approach, 
without a complete refit of all potential parameters. Obviously, the box subtraction will be obsolete in the future 
by the construction of a realistic potential which incorporates from the beginning explicitly nucleon and A d.o.f. in 
a coupled channel approach. Finally, we would like to remark that the inclusion of ( |88| ) in the Lippmann-Schwinger 
equation (pin leads to an energy dependence of the subtracted box graphs which we avoid by the substitution in (|3l|) 



R(z)V^ ] R(z) 



{R(z)V^ ] R(z)}\ z 



(90) 



F. The nucleon-nucleon interaction 



The properties of the two-nucleon system is governed by the NN interaction. In this work, we restrict ourselves 
to realistic NN potentials known from the literature, which have to be renormalized according to the discussion in 
the foregoing subsection. The starting points of our considerations are the Lippmann-Schwinger equation for the NN 
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scattering matrix ( |3l| ) and the Schrodinger equation for the effective, purely nucleonic component of the deuteron 
(|54|), which we write in the forms, respectively, 

T c Mz) = Vgf(z) + V^(z)G (z)T^(z) , (91) 
H + V^(M d ))\d)=M d \d), (92) 

with the effective A^A" interaction 

V$(z) = R(z) (V|L + [V Rx G (z)V xN ] con ) R(z) . (93) 

Note that in ( |9l"| ) the A isobar and the nd channel has been neglected, because the realistic NN potential models 
do not consider these additional degrees of freedom in general. Exceptions are, for example, the "full" Bonn and the 
Argonne V28 potentials |24[|33|| which include the A. 

In view of the discussion in the foregoing subsection, we follow different strategies in static and retarded calculations. 
In the static limit the effective interaction simplifies to 



for which we use the Bonn-OBEPR potential V OBEPR [|| to be renormalized according to Eq. (pi]), whereas in the 
retarded approach VSL in V^(z) has to be chosen as (cf. (|38|)): 

V^ = {Z$)^V*1\Z$)-1 (95) 
= - V Na (z)G (z)V A x(z)\ z=2Mn - V Nx G (z)T x (z)G (z)V xf j\ z=2MN , 

with 

Vn± = Z$ ] v££ m (Z$)- 1 (96) 



(note the analogy to ( p2|)) . Concerning the one-boson exchange part in (93), we are able to use the potential model 
yEister Q £ Ulster et al. ||20|pl[| , which can be considered as an extension of the retarded Bonn-OBEPT potential |24|] 
with respect to the inclusion of additional nN loops. Consequently, we use in the retarded calculations (note (BOh) 



V^(z) = V Els ^(z) + R^z)V^(z) , (97) 

Z- — 2X1 x 

where is given by (|9^), and the pionic part R K (z) of the dressing factor is discussed below. 

Because of its importance for the present work, it is worthwhile to study V Elster in some more detail. Concerning 
the one-boson exchange, all six mesons 7r, 77, a, 5, uj and p are taken into account whereas only the pion is considered 
in the dressing operator, therefore denoted by R 7r (z), for simplicity. This approximation is, however, not crucial with 
respect to unitarity up to the r\ threshold. Consequently, y Elster ( z ) has the structure 

V E * t * r (z) = fr(z) Y, \VN*G Q {z)V xR ] con R*{z) : (98) 

a;£:{7r,?7,er,($,ci;,p} 



where R 7r (z) is given in analogy to (33) with only irN loops in ZlL(z). Explicitly, one obtains for the matrix element 
of the latter 



with 



(NN; p'\Z( 2] (z)\NN; p) = S (p' - p) Zf 2] (z,p) (99) 



/ \2 (a ) 2 3 r°° k 2 r 1 

{zUz,p)) =l + 2 Ky ;> - - / dk—- dco S eF^((p-k) 2 ) 

\ W y ,y 'J 4tt Ane N {p) J e N (k) J_ x 

e N (p)e N (k) - M% -pkcosO 



Jn(jp- k) (e N (p) - e N (k) - tJ w (\p- k \j) (z - e N (k) - u>„(\jp— k \)j 
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where the angle between p and k is denoted by 9 

^ x by the o^v-i-c^uwi. 



In view of the renormalization of VR „ by the operator ZTL , renormalized meson- nucleon coupling constants 



fx 9° 
h = Z? 2] (E™,p) and 9x = Zfa(E™,p) (101) 

will appear in (|98l) , whose momentum dependence is small and, therefore, can be neglected (see Fig. |lO| ). The free 
parameters in (|9q), i.e., the cutoffs and physical coupling constants are fixed by fitting NN scattering data below it 
threshold and deuteron properties. The resulting values are presented in Tabic [jj 



IV. RESULTS FOR iViV SCATTERING 
A. Determination of the parameters 



Now we will fix th e remaining free parameters in the hadronic interaction by considering NN scattering. As dis- 
cussed in Sect. QIC, in a retarded approach the parameters of the diagr ams (a) through (c) of Fig. |q are in pr incipl e 



fixed by the parametrization of the meson-nucleon vertices (Sect. [II A) and of the ttNA vertex V^ A (Sect. Ill] 
Consequently, only the cutoffs and coupling constants in the diagrams (d') through (f ') can be treated as free param- 
eters. According to ( fFT| ) through (|80|), these diagrams contain five open parameters, namely the coupling constants 
f%N-n an d /aa?p' the parameter otANp and the form factors Fan Tr{q 2 ) and FANp(q 2 ) which are parametrized as in 
(|60|). However, the question arises whether the parametrization of the 7rAA vertex as obtained from fitting nN scat- 
tering data should be used for the OBE mechanisms, too. For example, one of the differences between NN and nN 
scattering is the fact that in NN scattering above n threshold a pion can be onshell, whereas in ttN scattering it must 
be onshell. This leads to dramatic differences in the cutoff Aajvtt in ((6^)- For example, the value of 1200 MeV (with 
"-A?r = 1) for AajVtt in the full Bonn potential |24| is much stronger than the value of typically 300 MeV (uatt = I ) 
or 500 MeV (uatt = 2) obtained from 7TA" scattering. Similar differences do occur also for the parametrization of the 
ttN vertex in NN versus irN scattering. 

Holzwarth and Machleidt |3^] have shown in a detailed analysis that these problems can be traced back to the 
use of the monopole or dipole parametrization of the hadronic form factors, which ties together the low and high 
momentum behaviour and thus is not able to describe both NN and nN data simultaneously with the same cutoff. 
In AW scattering, for example, pion momenta in the neighbourhood of q = are the relevant ones. Thus for small q, 
the form factor should be close to unity in order to achieve, for example, a quantitative description of the deuteron 
quadrupole moment, which is dominated by long-range mechanisms. Therefore, a large cutoff mass is needed. On the 
other hand, for large cutoff masses the form factor does not decrease fast enough with increasing pion momenta as is 
necessary in irN scattering. However, a form factor which is derived from the Skyrme model is able to describe both 
simultaneously, NN as well as irN scattering p6f , because it combines the features of a hard monopole or dipole form 
factor at q = with the quality of a soft form factor for larger pion momenta. 

For these reasons - although we are aware of the fact that this procedure leads to additional inconsistencies -, we do 
not use the values ( |67| ) in the OBE mechanism of our coupled NN-NA system. Thus, we replace the coupling constant 
f AN7T and the form factor Fantt{<1 2 ) in the retarded diagrams (a) and (b) of Fig. _|6| (see (|69| )) by the corresponding 
quantities fANTr($ 2 ) anc ^ -FAA»v(<f 2 ) of diagram (d') and (e'), respectively. For the A^A «-> A^A transition (diagram (c) 
and (f) in Fig. ||) we use the values ( |67| ) or ( |68| ) as obtained from ttN scattering. Similar to previous studies 
we do not determine the free parameters by a global fit to all NN channels in the A region. Instead, we concentrate 
on the most important 1 i?2 channel, because it is the only partial wave which couples to a NA-S state (sec Table ||) 
and thus where one expects the strongest NA interaction effects. 

In order to distinguish the various cases, we introduce as nomenclature "CC(approach, mesons)", where for "ap- 
proach" we consider a retarded one ("ret") and different static ones ("stat", "statl", and "stat2"). The entry "meson" 
refers to the mesons included in the A^A interaction, i.e., "7r" means only pion exchange and "n, p, oantt" means 
7r and p exchange including the additional coupling ( |76| ) whose strength is controlled by ctANTr-_ We begin with the 
discussion of the parameter choices for the retarded approach as listed in Table [II. For f AN ^, /aat p j and the cutoff 
Aaattt, we have chosen the values of the full Bonn Potential, i.e., 

(70 \ 2 (7o \ 2 

SlApLL = 0.224, KANp ' = 20.45, Aaa^*?* 2 ) = 1200 MeV, h A7T = 1 . (102) 
4-7T Ait 

It turned out that for a variety of different combinations of the remaining parameters a A Np and Aaa/> a satisfactory 
description of the 1 Z?2 phase shift is possible. Thus we have restricted the choice of otANp to the values —1, 0, 1. 
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The resulting cutoff masses Aaat p are listed in Table [II as determined by a fit to the 1 D 2 phase shift at Ti a j, rj 500 
MeV. The value a^Np = 1 corresponds to the usual neglect of the additional coupling (f76|) . On the other hand, the 
choice aANp = can be motivated by vector dominance insofar as for this choice the pNA vertex has the same spin 
structure as the dominant Ml 7 A" A coupling. 

As next, we turn to the static case. We would like to remind the reader that in this case all diagrams except the 
contribution (c) of Fig. ^ are incorporated into the static part of V° ^ , where all corresponding coupling constants (/°, 
fpi 9% f&Nin /ajv>)> the parameter a^Np and all form factors F„, Fantt and F&Np are treated as free parameters. 



The resulting values are listed in Table |VJ. In the simplest approach, CC(statl, 7r), which is identical to the one 
used in ]l3| , we neglect the ird channel and the p exchange in and completely, e.g., the A^A interaction is 

solely given by tt exchange, where the cutoff A^nn = A-ANn is fitted to the 1 D 2 partial wave at Ti a \, w 500 MeV. The 
ttNA vertex is given by the parametrization ( |68| ) of Sauer et al. ]2(|. Similar to the retarded case, we use these 
parameters ( |68[ ) also in the potential V^a (diagram (f) of Fig. ^). 

The second approach CC(stat2, tt) differs from CC(statl, tt) only with respect to the parametrization ( |67j ) instead 
of (^8|) for Va-, e.g., those values for are taken into account which are also present in our retarded approach. 
The choice ( p7| ) is also used in our third case CC(stat, tt, p, 0). But in contrast to the previous cases, we incorporate 

in addition p exchange in V^^} and V^A as wen as * ne K ^ channel. Furthermore, we set a/\Np = 0. 

B. Results 

We will start the discussion with the 1 D 2 partial wave. As already mentioned, all three retarded potential models 
are able to describe its phase shift and inelasticity equally well in the A region (see Fig. |ll|). This is of course not 
very surprising because we had used this channel for fixing the free parameters. However, above Ti a \, > 800 MeV 
one notes a rather large discrepancy between theory and experiment. In view of the fact that the parameter cxan p 
could not be fixed uniquely, we have studied the influence of different choices for a&Np on the other partial waves. It 
turned out that the strongest dependence was found for the P waves (see Fig. [j^). However, it was not possible to 
determine an optimal value for a^Np- While the 3 Po inelasticity and the 3 Pi channel seem to favour a^Np = 1> the 
value ctAjVp = — 1 leads to a slightly better fit of the 3 P 2 channel, especially for energies T; a f, < 700 MeV. But the 
overall description is still quite poor. On the other hand, this rather unsatisfactory situation is not very surprising 
since we have not fitted our parameters to all partial waves simultaneously. 

Now we will turn to the discussion of the various static and retarded approaches for the hadronic interaction. We 
have chosen the model CC(ret, tt, p, 0) as a starting point and compare it with the static approaches CC(statl, tt), 
CC(stat2, tt), and CCfstat, tt, p, 0) for all (t = 1) partial waves with total angular momentum j < 3 as shown 
in Figs. [13] through [l5|. With respect to the X D 2 phase shift, all four models yield more or less similar results. 
The inelasticity in CC(statl, tt) and CC(stat2, tt), however, is somewhat smaller than in the other two models. In 
accordance with [Q, this behaviour can be traced back to the neglect of the ird channel in CC(statl, tt) and CC(stat2, 
tt). For the other partial waves, the comparison with experimental data does not give a uniform picture. For example, 
the calculation CC(stat, tt, p, 0) yields the comparably best, but still not satisfactory description of the 3 Pi phase 
shift and the 3 Po inelasticity. On the other hand, the theoretical 3 P 2 phase shift is largely at variance with the data, 
whereas CC(statl, tt) and CC(stat2, tt) result in a rather good description. But in the 1 >S'o and 3 i 7 3 channels the 
retarded calculation CC(ret, tt, p, 0) is favoured. 

In summary, one finds that our fit procedure leads to a satisfactory agreement with the experimental data for the 
1 D 2 channel, which is the most important one in the A resonance region. The overall description of the other channels 
is fairly well but needs further improvements. As a first step one would need to construct a potential model (static 
or retarded) in which all parameters are fitted simultaneously to all partial waves for energies up to about T; a (, = 1 
GeV. This quite involved task will be one of our future projects. 

Another remark is in order with respect to the question whether the static or the retarded interaction gives a better 
description. Although in principle, from the basic physical ideas, a retarded interaction is more appropriate, our 
results for the pure hadronic reactions do not allow a clear answer. To this end, one has to consider e.m. reactions 
where first results Ei 15 indicate a preference for the retarded interaction. This will be the subject of forthcoming 
papers. 

C. Comparison with other approaches 

In the coupled NN-NA approach of Lcidcmann et al. |l2| , the static Reid soft core and Argonne potentials have 
been used as starting point for the NN interaction, whereas the Trd channel has been neglected. Since the numerical 
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calculations have been performed in r-space, nonlocalities, which are, for example, present in the TVA propagation, 
could not be taken into account exactly. Therefore, additional approximations in the box-renormalization procedure 
had been used which lead to a much more flexible treatment by introducing one additional free parameter for each TVTV 
partial wave. Thus, compared with our results, a better description of TV TV scattering data was achieved because of 
this additional phenomenological ingredient. With respect to the static framework, we further would like to mention 
the work of Wilbois jj7|. Among other things, he compared the results of a coupled channel approach for different 
static NN potentials (Bonn-OBEPR, Bonn-OBEPQ, Paris and Nijmegen). It turned out that the TVTV phase shifts 
and inelasticities are rather sensitive to the choice of the underlying TV TV potential. No TV TV potential, incorporated 
into the CC approach, has been able to produce a quantitative description of all TVTV phases and inelasticities in the 
A resonance region. 

Retarded approaches based on three-body theory [ p7|p8| have already been discussed in Sect. |. Here, we want to 
consider in some more detail the work of Bulla et al. and Elster et al. ]2^ , |23| ] . The coupled channel approach of p3[ , 
which includes an explicit 7rTV vertex, differs at least in three important aspects compared to the present approach: 
(i) The ird channel is neglected, (ii) Retardation is treated only for pion exchange. For the cutoff mass of the corre- 
sponding 7rTV vertex form factor, a very small value of 443 MeV has been used in dipole parametrization as obtained 
from a fit of 7rTV scattering data in the Pu channel, whereas the value of 1700 MeV in the Elster potential results 
from a fit to TVTV scattering data, (iii) The static Paris potential is used as the underlying TVTV interaction. Therefore, 
because of the renormalization procedure (^9|) , the effective retarded OPE potential is considerably weakened. In view 
of the points (ii) and (iii), it is not very surprising that Bulla et al. found only a small influence from retardation, in 
fact much smaller than in our approach. Moreover, due to the mixture of retarded and static frameworks it is rather 
questionable to apply this approach to electromagnetic reactions on the deuteron. For example, due to the existence 
of two different tt exchange mechanisms in the Paris potential and the retarded interaction, two different types of 
MECs have to be constructed in order to guarantee gauge invariance - which is rather artificial. 

In the work of Elster et al. |22j , who have extended the Bonn boson exchange potential into the region above pion 
threshold by the inclusion of the A isobar, the interactions are treated in a retarded approach including appropriate 
nucleon and A dressing. Elster et al. start within a Lee model framework whereas our retarded interactions are 
derived within time-ordered perturbation theory. In contrast to our model, the ird channel as well as the TVA <-> TVA 
transition potential (diagram (c) and (f ) in Fig. @) are neglected, so that unitarity is violated above 7r threshold. On 
the other hand, the A is treated relativistically whereas we describe it nonrelativistically. Moreover, AA components 
are taken into account in [p2| . Concerning the pure nucleonic sector, however, both approaches yield exactly the same 
interaction. 



V. SUMMARY AND OUTLOOK 



In this paper, we have presented a hadronic interaction model which is suited for the study of electromagnetic and 
hadronic reactions in the two-nucleon sector for exitation energies up to about 500 MeV, for example, TVTV scattering 
or electromagnetic deuteron break-up, in which not more than one pion is created or absorbed. This model respects 
in particular three-body unitarity. It is based on a previously developed model of Sauer and collaborators |p3| , p5| , p6| 
which allows the explicit consideration of retardation in the two-body meson-exchange operators within a TVTV-TVA 
coupled channel approach using time-ordered perturbation theory. 

Since retarded interactions are not hermitean, we have generated the retarded one-boson exchange mechanisms by 
considering explicitly meson-nucleon and 7rTVA vertices. Therefore, additional mesonic degrees of freedom besides the 
baryonic components had to be included explicitly into the Hilbert space of a two-baryon system. In the present work, 
we have restricted ourselves to the one-meson-approximation, which means that we allow only configurations with 
one meson present besides the baryons. As mesons we have taken into account explicitly n, p, u>, a, 5, and r\. In order 
to satisfy two- and three-body unitarity, we have incorporated the ~nd channel as well as intermediate pion-nucleon 
loops. In view of the latter mechanism, a distinction between bare and physical nucleons was necessary in order to 
avoid inconsistencies. 

In our explicit realization within a field-theoretical framework, we have used the realistic retarded potential model of 
Elster et al. | 20|pl"[ ] as input for the basic nucleon-nucleon interaction in pure nucleonic space. This potential had to be 
renormalized because of the additional d.o.f., which lead to further contributions from TVA and ird states. Because of 
some necessary approximations in our numerical evaluations, unitarity is not completely obeyed, which fact, however, 
we do not consider very critical. The free parameters of our model have been fixed by fitting 7rTV scattering in the P33 
channel and TVTV scattering in the 1 £>2 channel. For practical reasons, no global fit of the parameters to all relevant 
TVTV or 7rTV channels has been performed. Therefore, the overall description ot the TVTV phase shifts and inelasticities 
is fairly well but needs some further improvement in the future. Moreover, due to our fit procedure, there remain 
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some ambiguities because not all parameters could be fitted uniquely. For example, besides the usual pNA interaction 
density ([75]) we have incorporated an additional alternative ([76|) into the p exchange of the NN — > NA transition, 
which has, at least to our knowledge, not yet been discussed in the literature. The corresponding coupling constant 
could not be fixed unambiguously by studying the partial wave alone. However, different choices of this coupling 
constant lead to rather drastic changes in the P wave phase shifts in the A region. For an improved description one 
would need to construct from scratch a hadronic interaction model which incorporates from the beginning nucleon, 
meson, and A degrees of freedom and whose open parameters are fitted to the phase shifts and inelasticities of all 
relevant NN scattering partial waves for Ti a b energies up to about 1 GeV. 

Another considerable improvement would be the implementation of the convolution approach of Kvinikhidze and 
Blankleider |4C|-|4^] into our model. This would allow one to abandon the one-meson-approximation which causes 
several pathologies. From a conceptual point of view, this problem deserves further attention because it deals with 
the fundamental question of deriving off-shell properties of particles in a many-body system (for example the dressing) 
from those of a single-particle system without creating inconsistencies. 

In conclusion, we believe that the present model is realistic enough for the study of electromagnetic reactions like 
photo- and electrodisintegration of the deuteron as well as pion production. The predictions for the observables of 
such reactions will be presented in forthcoming papers. 



[1] C. Herberg et al., Eur. Phys. Journal A 5, 131 (1999); M. Ostrick et al., Phys. Rev. Lett. 83, 276 (1999). 
[2] I. Passchier et al, Phys. Rev. Lett. 82, 4988 (1999). 

[3] H. Arenhovel, G. Krefi, R. Schmidt, and P. Wilhelm, Phys. Lett. B 407, 1 (1997); Nucl. Phys. A 631, 612c (1998). 
[4] R. Machleidt, Adv. Nucl. Phys. Vol. 19, p. 189 (Plenum, New York, 1989) 
[5] J. Carlson and R. Schiavilla, Rev. Mod. Phys. 70, 743 (1998). 
[6] H. Arenhovel and M. Sanzone, Few-Body Syst. Supp. 3, 1 (1991). 
[7] H. Arenhovel, Few-Body Syst. 26, 43 (1999). 

[8] H. Arenhovel, F. Ritz, and Th. Wilbois, Phys. Rev. C 61, 034002 (2000). 

[9] J.M. Laget, Nucl. Phys. A 312, 265 (1978); Can. J. Phys. 62, 1046 (1984). 
[10] M. Schwamb, H. Arenhovel, and P. Wilhelm, Few-Body Syst. 19, 121 (1995). 
[11] H. Tanabe and K. Ohta, Phys. Rev. C 40, 1905 (1989). 
[12] W. Leidemann and H. Arenhovel, Nucl. Phys. A 465, 573 (1987). 
[13] P. Wilhelm and H. Arenhovel, Phys. Lett. B 318, 410 (1993). 

[14] M. Schwamb, H. Arenhovel, P. Wilhelm, and Th. Wilbois, Phys. Lett. B 420 , 255 (1998). 
[15] M. Schwamb, H. Arenhovel, P. Wilhelm, and Th. Wilbois, Nucl. Phys. A 631, 583c (1998). 
[16] M. Schwamb, PhD Thesis, Mainz 1999. 

[17] W.M. Kloet, R.R. Silbar, Nucl. Phys. A 338, 281 (1980); R.R. Silbar, W.M. Kloet, Nucl. Phys. A 338, 317 (1980). 

[18] H. Tanabe, K. Ohta, Nucl. Phys. A 484, 493 (1988). 

[19] W.M. Kloet, R.R. Silbar, Phys. Rev. Lett. 45, 970 (1980). 

[20] Ch. Elster, PhD Thesis, Bonn 1986. 

[21] Ch. Elster, W. Ferchlander, K. Holinde, D. Schiitte, and R. Machleidt, Phys. Rev. C 37, 1647 (1988). 

[22] Ch. Elster, K. Holinde, D. Schiitte, and R. Machleidt, Phys. Rev. C 38, 1828 (1988). 

[23] A. Bulla and P.U. Sauer, Few-Body Syst. 12, 141 (1992). 

[24] R. Machleidt, K. Holinde, and Ch. Elster, Phys. Rep. 149, 1 (1987). 

[25] P.U. Sauer, Progress in Particle and Nuclear Physics, Vol. 16 , ed. A. Faessler (Pergamon Press, 1986). 
[26] H. Popping, P.U. Sauer, and X.-Z. Zang, Nucl. Phys. A 474, 557 (1987); 

H. Popping, P.U. Sauer, and X.-Z. Zang, Nucl. Phys. A 550, 563 (1992). 
[27] M. T. Pena, H. Garcilazo, U. Oelfke, and P. U. Sauer, Phys. Rev. C 45, 1487 (1992). 
[28] P. Wilhelm and H. Arenhovel, Nucl. Phys. A 593, 435 (1995). 
[29] P. Wilhelm and H. Arenhovel, Nucl. Phys. A 609, 469 (1996). 

[30] M. Lacombe, B. Loiseau, J. M. Richard, R. Vinh Mau, J. Cote, P. Pires, and R. de Tourreil, Phys. Rev. C 21, 861 (1980). 
[31] R.A. Arndt et al, program SAID. 

[32] H. Goller and H. Arenhovel, Few-Body Syst. 13, 117 (1992). 

[33] R.B. Wiringa, R.A. Smith, and T. L. Ainsworth, Phys. Rev. C 29, 1207 (1984). 

[34] V.G.J. Stoks, R.A.M. Klomp, C.P.F. Terheggen, and J.J. de Swart, Phys. Rev. C 49, 2950 (1993). 

[35] A.M. Green and M.E. Sainio, J. Phys. G. 8, 1337 (1982). 

[36] G. Holzwarth and R. Machleidt, Phys. Rev. C 5, 1088 (1997). 



17 



[37] Th. Wilbois, PhD Thesis, Mainz 1996. 

[38] H. Tanabe and K. Ohta, Phys. Rev. C 36, 2495 (1987). 

[39] T.D. Lee, Phys. Rev. 95, 1329 (1954). 

[40] A.N. Kvinikhidze and B. Blankleider, Phys. Rev. C 48, 25 (1993). 

[41] A.N. Kvinikhidze and B. Blankleider, Phys. Lett. B 307, 7 (1993). 

[42] B. Blankleider and A.N. Kvinikhidze, Few-Body Syst. Suppl. 7, 294 (1994). 

[43] P.U. Sauer, M. Sawicki, and S. Furu, Prog. Theor. Phys. 74, 1290 (1985). 



APPENDIX A: THE PHYSICAL NUCLEON 



Here we will consider briefly a "physical" nucleon state \N). It is straightforward to show that the physical nucleon 
state can be written as 

\N(p)) = N^(p) (1 + go(e N (p))v° x s) \N(p)) , (Al) 
where the free propagator go(z) is defined by 

g (z) = {z- ho)- 1 , (A2) 

and v xR is the only interaction term describing the emission of a meson, since we do not consider a diagonal meson- 

nucleon interaction in , i.e., v xx vanishes. Note that a A will not appear because of isospin conservation and the 
one-meson-approximation. For the renormalization constant one finds 



N [1] (p) = y Jl + vf 1] (p), (A3) 

where (p) is defined by 

(N(f)\v NX go(e N (p'))go(e N (p))v X x\N(p)) = ^{p)5{f - p) . (A4) 
Furthermore, for the counter term one obtains the identity 

(N(p')\v%\N(p)) = -(N(P)K X (e N ( P ) - ^xx)-V^|iV(p)><5(i?' -p) , (A5) 
yielding the following compact expression 

dzS(z-H )v% x g (z)v xff . (A6) 

In summary, it is obvious that a distinction between bare and physical nucleons is necessary due to the occurrence of 
the interaction v% x which generates intermediate meson-nucleon loops. 



APPENDIX B: DERIVATION OF THE NN SCATTERING T-MATRIX 



In this appendix we will briefly sketch the derivation of Eq. (p0[). With the help of projection operators one can 
bring ( p7|) into the equivalent form 

\NN; p, a) (±) = N^(f) (l + G„« w ± ie)T pR {E XN ± ie) + G Q (E XN ± te)T xR (E XN ± ie)) \NN; p, a) , (Bl) 

where T pp and T xp are given by (note P = P R + Pa) 

T° PP (z) = Vj> p (z) + V^ P (z)G Q (z)T a PP (z) , (B2) 
T° XP (z) = G^(z)G x (z)V XP (l + G (z)T PP (z)) . (B3) 

In these relations, V PP is defined by 

V$ P (z) = V^ P + V^ x G x (z)V XP , (B4) 
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and G x {z) in (pCf). The driving term V PP in (B4) can now be split into a connected ("con") and a disconnected 
("dis") part 

V° P {z) = V°f s (z) + V°p Con (z) , (B5) 

where, by definition, V pp dls contains only those parts of the driving term V PP which do not describe an interaction 
between the two baryons except for the ttN loop contributions to the A self energy, which we hav e in corporated into 
Vp j3 CO "(.z), containing otherwise the genuine baryon-baryon interaction. In detail one finds from (B4) and (|5p| ) with 

V°f*(z) = V® R + [V% x Go(z)V^] dis , (B6) 
V°>™(z) = V°P + [V^ x G Q (z)V%] con + V^ x G (z)T x (z)G a (z)V% + [Vl x G (z)V XA ] ^ , (B7) 

where we have defined 

KxG (z)V Xa ] d . s = J2 <x(i)G (z)v Xa (i] , a £ {N, A} , (B8) 

i=l,2 

K x G Q {z)V x0 ] con = vZx(})Go(z)v%(j), a, P £ {N, A} . (B9) 

Thus the connected part V PP con (z), which is shown in Fig. ||, contains besides V PP ^ a retarded one-boson exchange 
interaction, the coupling to the nNN channel, and, furthermore, also a disc onnected part, namely, the already 
mentioned ttN loop contributions to the A self energy in the last term of (B7). The separation ( |B5| ) allows us to 
represent the total amplitude by a "disconnected" T pp is (z) and a "connected" T°' con (z) amplitude 

T° PP (z) = T° p f s (z) + (l + T° p f s (z)G (zj) T° P ™ n (z) (l + G (z)T° p f s (zj) , (BIO) 

where 

T° p f s (z) = V°f s {z) + V^ dls (z)G Q (z)T^ dls (z) , (Bll) 
T% on {z) = V^ on {z) + V°' p con (z)G (z)T^ on (z) , (B12) 

with the dressed propagator 

G (z) = (z-H -V°f s (z))- 1 

= G (z) + G (z)T° P f s (z)G (z) . (B13) 

A graphical representation of the dressed propagator in the form Gq(z) = Go(z) + Gq(z) V pp dls (z)Go(z) is displayed 
in Fig. [l6| and of the connected scattering amplitude T pp on (z) in Fig. 
We are now in the position to rewrite (E7f) as follows 



\NN; p, a)& = N^(p) (1 + G X (E™ ± ie)V XP ) (l + G (E™ ± i e)T° p r(E? N ± ie) 

x (l + Go(E XN ± ie)T° p f s (E™ ± ie)) \NN; ft a) . (B14) 



In order to determine the counterterm V^k in the two-nucleon sector, we will consider the special case of two 
noninteracting physical nucleons denoted by \NN; p, a). Using the decomposition 

\NN; p,a) = Pn\NN; p,a) + P x \NN; p,a) (B15) 

into nucleonic and mesonic components, one obtains, for example, for Pfj\NN] p, a) a Schrodinger equation 

(H + V^f s (E XN )) P„\NN; p, a) = E^ N \NN; p, a) . (B16) 

In view of Pn\NN) ~ \NN) and the trivial identity H P^\NN(p, a)) = E^ N P^\N N (p, a)) , one finds for the 
counterterm (note the analogy to the one-nucleon sector) 
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NN 



dzS(z-H ) [V$ x G (z)V° 



Nidi 



(B17) 



It should be emphasized that, due to the one-meson- approximation, one has to face the pathological situation that 
a dressing of both nucleons at the same time is not possible. Therefore, \NN; p, a) is not the direct product of two 
physical nucleon states of the one-nucleon sector, i.e., 



\NN(p)) ? \N(p)) ® \N(-p)) 



(B18) 



A way out of this well known problem EjJ] could be the "convolution approach" of Kvinikhidze and Blankleider 
E which for technical reasons has not yet been implemented in our model. 

With the help of the two-body renormalization operator Z[ 2 ](z) as defined in (|3~4|), one can rewrite the dressed 
propagator Gq(z) according to 



G (z) = G (z)Zpf(z), 



(B19) 



and one obtains from ( |B13| ) the identity 

l + G (z)T° p f s (z) = Z^(z), (B20) 
which allows us to determine the renormalization constant Nn] (p) straightforwardly. From the normalization condition 

(NN;p\a'\NN;p,a)=5 a , a 6(p' -p), (B21) 

one finds 

(NN; p',a\{Z$)- x \NN; p,a) = N [2] (p)6 a , a 5(p' - p) . (B22) 

Introducing the "renormalized" interaction from (B2|), one obtains from JB12| ) with the help of ( pig| ) and @, for 
the renormalized amplitude T p ° p (z), defined as 

T??{z) = %\{z) T% m {z) , Z^{z) (B23) 



the Lippmann-Schwinger equation in (|31j), which is very useful for practical evaluations. In fact, considering NN 
scattering for which the T matrix is given by the onshell matrix element of T° 



Npfo) (NN; f, a'\T°(E^ N + ie)\NN; p, a) , 



where \p'\ = \p \ = p has been used, one can rewrite the total amplitude Tp P (BIO) as 

T PP (z) = T° p f s (z) + Z- 2 \(z) T&(z) Z^(z) . 



(B24) 



(B25) 



With (|B22|) and the relation 

T° P ^ s (z)\NN; ft a) = (z^f{z) - l) G^(z)\NN; p,a) , (B26) 
which vanishes for z = E^ N ± ie, one finds for the matrix element in Eq. ( |B24 ) the desired expression of (pC|). 



TABLE I. Meson parameters of V Elster . The parameters for the a meson apply only to the (t — 1)-NN potential. For t = 0, 
we use gl/4n = 9.4050, m CT = 580 MeV, A CT = 2300 MeV, n a = 1. 



iL(/»/s. 



[MeV] 



A, [MeV] 



7T 

9 

U! 

a 
6 
V 



14.4 
0.9 (6.1) 
20.0 (0) 
9.4080 
0.3912 
5.0 



138.03 
769 

782.6 
575 
983 

548.8 



1700 
1500 
1500 
2600 
1500 
1500 
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TABLE II. NN and NA partial waves ( lj) with total angular momentum j < 3, parity n and isospin t. 



j 


7T 




t 




NN 




JVA 






+ 




1 
1 




L S 
3 Po 




b D 
3 Po 


1 
1 
1 


+ 




1 






3 Pi 
3 S U 3 D 1 
x Pi 




J Pi, °Pi, ;, Pi 


2 
2 
2 


+ 
+ 




1 
1 






L D 2 

3 P 2 , 3 F 2 
3 D 2 




°S 2 , *D 2 , a D 2 , a G 2 
3 P 2 , 5 P 2 , 3 F 2 , 5 F 2 


3 
3 

3 


+ 




1 







:i F 3 
3 D 3 , 3 G 3 




'/'-.. 'F;. •/•,. II., 




TABLE III. Parameter values for V$ A 


and Vaa hi the retarded approach. 




CC(type) 


(■^AJVtt) 
4-7T 


Aaat^ [MeV] 


riAw 


(fO \2 


Aaivp [MeV] 


flAp CtANp 


(ret, 7r, p, —1) 
(ret, 7r, p, 0) 
(ret,7r,p, 1) 


0.224 
0.224 
0.224 


1200 
1200 
1200 




1 
1 
1 


20.45 
20.45 
20.45 


1340 
1420 
1560 


2 -1 
2 
2 1 


TABLE IV. Parameter values in 
Vatt ( see the discussion in the text) 


the static approach. Note that the models differ 


in the parametrization of the tvNA vertex 


CC(type) 


Hi) 2 

4-7T 


K [MeV] 






4tt {fp ! 9p) 


A p [MeV] 


rip 


(statl, 7r) 
(stat2, 7r) 
(stat, 7r, p, 0) 


0.08 
0.08 
0.0778 


700 
680 
1300 




1 
1 
1 


0(0) 
0(0) 
0.84 (6.1) 


1400 


1 


CC(type) 




A A iv,r [MeV] 






{f ANp) 


A ANp [MeV] 


n Ap CtANp 


(statl, 7r) 
(stat2, 7r) 
(stat, 7r, p, 0) 


0.35 
0.35 
0.224 


700 
680 
1200 




1 
1 
1 


20.45 


1600 


2 
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FIG. 1. Diagrammatic representation of the various components of V° . The open ellipse symbolizes a given hermitean 
two-body interaction. The one-nucleon counter term v [c] is indicated by a cross. 
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FIG. 2. Diagrammatic representation of the free propagator Go(z) with the three contributions in TL^ , Ti A , and Ti,\ 
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FIG. 3. Diagrammatic representation of the connected driving term Vp P con (z). The greek letters a and a' label either a bare 
nucleon N or a A. 
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FIG. 4. Diagrammatic representation of the iViV scattering matrix T x (z) in the presence of a spectator meson. 



FIG. 5. Diagrammatic representation of the dressed A propagator g& (see Eq. (|53|)). In order to distinguish g& from the 
free propagator, it is denoted by a 
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(a) (b) (c) 



(d) (e) (f) (d') (e') (f) 

FIG. 6. One-boson exchange diagrams for the various interactions involving a NA configuration. Diagrams (a) through (c) 
represent retarded contributions (here it only) . Diagrams (d) through (f ) are treated in the energy independent approximation 
as indicated by the diagrams (d') through (f ') because of the absence of 7riVA and 7rAA configurations. 




(a) (b) 

FIG. 7. Diagrammatic representation of (a) V X x describing the interaction V NN between the two nucleons in the 3 'S\- Z D\ 
channel with a pion as spectator, and (b) of the term [Vpx Go(z) T x (z) Go(z) Vxp]con in Eq. (fi(j|). 
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FIG. 8. Phase shifts and mixing angle for the 3 Si- 3 D\ channel of NN scattering as a function of the lab kinetic energy Ti a i,, 
compared with the experimental data (solution SM97 of Arndt et al. Notation of the curves: dotted: separable ansatz 

(|^) with the deuteron wave function of the Bonn OBEPR potential; solid: exact calculation using the full OBEPR potential. 
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(a) 



(b) 



FIG. 9. Dispersive contributions to the AW interaction: (a) from intermediate iVA states, (b) from intermediate offshell 
NN scattering in the presence of a spectator meson. 
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FIG. 10. The renormalization factor (Ep N , p) of Eq. (JQS|) as function of the external nucleon momentum p for different 



(s°Y (a ) 2 
parameter sets: dashed: v 7' — 15, = 500 MeV; dash-dotted: v 2 — 15, A„ 



1700 MeV; dotted: 
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MeV; solid: 



W) 2 



= 25, A, 



1700 MeV. A dipole form factor has been used. 
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FIG. 11. Phase shift 5 and inelasticity p for the 1 I?2 channel in comparison with experiment (solution SM97 of Arndt et 
al. | ^H ) for different retarded potential models (see Table III): full curve: CC(ret, ty, p, 0), dotted curve: CC(ret, 7T, p, 1), and 
dashed curve: CC(ret, n, p, — 1). 
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FIG. 12. Phase shift 5 and inelasticity p for the 3 Po, 3 Pi, and 3 P2 channels for different retarded potential models as in 
Fig. [y]. Notation of the curves as in Fig. [y] and experimental data: solution SM97 of Arndt et al. (3lj . 
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FIG. 13. Phase shift S and inelasticity p for the 1 S'o , 3 Po and 3 Pi channels for different static and retarded approaches: 
dotted curve: CC(statl, n), dashed curve: CC(stat2,7r), dash-dotted curve: CC(stat, n, p, 0), and full curve: CC(ret, n, p, 0). 
The experimental data represent solution SM97 of Arndt et al. |3l|. 
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FIG. 14. Phase shift S and inelasticity p for the x Di and 3 f3 channels for different static and retarded approaches as in 
Fig. [13. Notation of the curves as in Fig. [13 and experimental data: solution SM97 of Arndt et al. pl[ . 
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FIG. 15. Phase shifts, inelasticities and mixing angles for the 3 P2- 3 ' F2 channel for different static and retarded approaches 
as in Fig. [l^. Notation of the curves as in Fig. [lj] and experimental data: solution SM97 of Arndt et al. J31] . 
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FIG. 16. Diagrammatic representation of the dressed propagator Go{z). The cross represents the one-nucleon counter term 
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FIG. 17. Graphical representation of the scattering amplitude Tpp 0n (z). The greek letters a and a' label either a bare 
nucleon TV or a A. The driving term Vp'p° n (z) is shown in Fig. |J| 
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